Interactions between colloidal particles induced by polymer brushes grafted onto the 

substrate 
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We investigate the interaction energy between two colloidal particles on or immersed in nonad- 
sorbing polymer brushes grafted onto the substrate as a function of the separation of the particles 
by use of self-consistent field theory calculation. Depending on the colloidal size and the pene- 
tration depth, we demonstrate an existence of repulsive energy barrier of several ksT, which can 
be interpreted by separating the interaction energy into three parts: colloids-polymer interfacial 
energy, entropic contribution due to "depletion zone" overlap of colloidal particles, and entropically 
elastic energy of grafted chains by compression of particles. The existence of repulsive barrier which 
is of entirely entropic origin, can lead to kinetic stabilization of the mixture rather than depletion 
flocculation or phase separation. Therefore, the present result may suggest an approach to control 
the self-assembling behavior of colloids for the formation of target structures, by tuning the colloidal 
interaction on the grafting substrate under appropriate selection of colloidal size, effective gravity 
(influencing the penetration depth), and brush coverage density. 



Assembling colloidal nanoparticles into ordered struc- 
tures is essential to the fabrication of nano-materials with 
high performance in optical, electrical, and mechanical 
properties. [1111110 iJ,[ia[Il[ll[ll[li 
Due to the small length scale and interaction energy, the 
ordered structures normally are realized by self-assembly 
of nanoparticles in solution or at an interface for the re- 
quired motion to reach minimum energy state. The 
mechanism of such self-assembling behavior greatly relies 
on the effective interaction between particles mediated by 
the environment such as solvent, substrate, and external 
field. The fundamental element is obviously the effective 
interaction between two individual colloidal nanoparti- 
cles as a function of their distance, and the results will 
be helpful to our understanding of the collective assem- 
bling behavior of large amount of particles. Usually, the 
interaction between large colloidal particles in a solution 
of small ones or nonadsorbing polymers is strongly influ- 
enced by the entropic excluded- volume (depletion) forces 
which were first recognized theoretically by Askura and 
Oosawa(AO) . (lfij Such an effect can be successfully used 
to explain phase separation in colloid-colloid and colloid- 
polymer mixtures, [ljj Interestingly, Lin et al.[13j com- 
bined entropic depletion and patterned surfaces to study 
the self-assembly of colloidal spheres on periodically pat- 
terned templates. Many theoretical and experimental 
works have already been performed to examine the en- 
tropic depletion interactions between colloidal particles 
in different suspensions as depletion agents. The evolved 
systems include colloid/polymer mixture, fTA fTst IT^. |20| 
colloids in DNA solution s l21l o r nematic liquid crys- 
tal/rigid rod backgrounchjpajp, H HI |H HI [H con- 
fined colloidal systems |29l l3d|| and binary hard-sphere 
mixtures. |3ll 1321 |33j In colloid/polymer suspensions, 
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when the distance between sphere surfaces or between 
sphere surface and the wall is smaller than the diame- 
ter of polymer coils, there will be an attractive potential 
which depends strongly on the polymer concentration. 
The topology of phase diagram depends on the ratio of 
polymer radius of gyration to the radius of colloidal par- 
ticles. Usually, in the colloid limit with small polymer 
to colloid size ratios, depletion interactions can be well 
described in pairwise form. On the contrary, the pro- 
tein limit for larger polymer sizes requires the incorpo- 
ration of many-body contributions. [20] Direct measure- 
ment of depletion of colloids in semidilute DNA solutions 
has revealed the range of "depletion zones" proportional 
to the correlation length. |2l| Recently, Fredrickson et al. 
|34j investigated the potential of mean force between two 
nanoparticles in a symmetric diblock copolymer matrix 
using self-consistent-field (SCF) theory. 

On the other hand, the grafting of polymer chains to 
surfaces is a widely used method to modify their proper- 
ties including adhesion, lubrication, and wetting behav- 
ior, and has many useful applications such as colloidal 
stabilization, polymeric surfactant, biocompatibility, and 
drug carriers. SCF theory formulated in bispherical co- 
ordinates was developed to study the interaction between 
surface-grafted spherical assemblies, and the results show 
hints to the conditions for colloidal stabilization. [35l |3(J 
Recently, polymer-coated surface has become an alterna- 
tive substrate for the assembly of colloidal nanoparticles 
and metal nanocrystals because of the "softness" prop- 
erty of interfaces and the nanoscale characteristic length. 
The entropic elastic energy of such polymer brushes will 
play an important role in colloidal assembly since the 
typical energies of both polymer-chain stretching and col- 
loidal assembly are comparable to thermal energy. For 
example, the grafted polymer always exerts a repulsive 
entropic force on incoming particles, and most works fo- 
cused on the interaction between polymer brushes and 
individual incoming particles and how to prevent the ad- 
sorption of colloids such as proteins onto surfaces under 
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FIG. 1: A schematic of the two colloidal particles immersed 
in polymer brushes. 



various grafting density, chain length, and interactions 
between chains and surface. |37l IM IM lit! lil lil 
Some experiments have adopted polymer-grafted surface 
to control the assembly of nanoparticles. |E 0> H> El As- 
semblies of nanoparticles with number density gradients 
in two and three dimensions are achieved. However, to 
the best of our knowledge, there have been no systematic 
theoretical and experimental studies into the interaction 
energy between two colloidal particles on or immersed in 
the polymer brushes, which is vitally important to un- 
derstanding the assembling behavior of many colloidal 
particles on the polymer-grafted substrate. 

In the present paper, we propose a theoretical model to 
study the effective interaction between two hard colloidal 
particles on or immersed in polymer brushes with differ- 
ent penetrating depths and particle's radii comparable 
to the natural size of polymer coils. We try to present 
a detailed and complete picture on such unique interac- 
tion potential, in order to understand the kinetic stabi- 
lization^. e., colloidal dispersions) and flocculation of col- 
loids due to the presence of grafting polymer chains. The 
results are expected to be helpful to assembling experi- 
ments of particles on such a substrate with end-grafted 
polymers. 

In our two-dimensional model, the system is placed 
in xz plane and the substrate is horizontally fixed at 
z = with n a uniformly grafted polymer chains (see 
Fig.l). On the top of the brush is the polymer sol- 
vent composed of free chains. Two colloidal particles 
are modelled by hard circles of radius R located a dis- 
tance d + 2R apart and a distance h from the substrate; 
|44j there are no attractive or repulsive interactions be- 
tween them. Each colloidal particle is surrounded by a 
depletion zone of thickness A, where the concentration 
of polymers becomes depleted compared to that of the 
bulk solution. ^3,0] All polymer chains and colloids are 
allowed in the region < z < Z max . Here, Z max is large 
enough to avoid the influence of brushes and particles, 
and as Z approaches Z max , the free "solvent" chains at- 
tain their bulk properties. 45] The volume of the system 
V is L x x Z max , where L x is the lateral length of the 
surface along the x axes. The grafting density is defined 
as a = n a / L x . All polymer chains are of the same poly- 
merization index N and flexible with the same statistical 
length a, and incompressible with a segment volume p^ 1 . 



The probability distribution for molecular conformations 
of a Gaussian chain a is assumed to take the Wiener form 
P[r a (s)] cc exp[-^j f Q ds\ ^ff^- \ 2 ], where r a (s) denotes 
the position of segment s on chain a. In every calcula- 
tion, the positions of two colloids are fixed, and the local 
volume fraction <p p (r) of the colloidal particle of radius 
R with center at position r c is described by taking into 
account the depletion interactional! 0] between colloids 
and polymer chains, 
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|r-r c | <R, 
R < |r - r c | < R -| 
R + A < |r - r r 



We take the incompressibility constraint that y a (r) + 
ipp (r) + (fip(r) — 1.0, which will ensure the exclusion of 
polymer chains from the depletion zones surrounding the 
two colloidal particles. Here, p tt (r) and <pp{v) repre- 
sent local volume fractions of grafted and free chains, 
respectively. The size A of depletion zone is deter- 
mined by depletion agents in the surrounding environ- 
ment, and is independent of the radius R of colloids. 
For highly concentrated polymer solutions the "deple- 
tion zone" thickness is chosen to be about the size of a 
free segment plj. and here we set A = a.fis) We define 
Flory-Huggins interaction parameters between brush- 
solvent(xa/3), brush-particle {xap), and solvent-particle 
(X0 P ), and take XapN = 0.0 and Xa P N = Xf) P N = 20.0 
since we assume that all the polymers have the same 
chemical nature, while the colloidal particles are insol- 
uble to polymers. We use the grand canonical form 
of SCF theory, which has been proven to be power- 
ful in calculating equilibrium morphologies in po l yme ric 

system, |H E ElUmiHIUHlEllH ^ ^ t0 

deal with polymer solvents and brushes. 

The grand canonical partition function for the sys- 
tem can be written as 



e nf,Nf *Z„ a 



(1) 



n B =0 



where /z is the chemical potential per segment of solvent 
chains. Z nottn „ is the canonical partition function for n a 
grafted chains and np free polymers: 



Zn a ,n, = —r^ / J] Dr<*(s)P[r a (s)] n Dr (s) 



/3=1 



P[ra(s)]<5[l -<p a — tpp- ifip] exp[ 
lp(r?(0)) > 

a=l 



V 



k B T' 



(2) 



where the first S function enforces incompressibility and 
the second ensures the anchoring of the chain ends on 
the substrate, kg is the Boltzmann constant, T is the 
temperature, and v is the interaction energy. (p a and 
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(pp are corresponding operators to local concentration ip a 
and ipp. The end of brush chains can move on the sub- 
strates, although the total number of chains on surfaces is 
fixed(we take a — 0.25). This is so-called liquid brushes, 
contrary to solid brushes where the immobile chains are 
anchored onto surfaces^^- Here we consider the liquid 
brush case because of its wide-ranging applications in col- 
loidal self-assembly and biological organization including 
cell adhesion and interactions between polymer-coated 
membranes and proteins(or cells) |4jj. The SCF theory 
gives the free energy 



iVF 



p k B TV 



<pV p V 



Np, 



Nv 



p a k B TV 

-l/V J dr[£(l -tp a -<pp- P P ) 

+w a ip a + wpepp] , (3) 

where 4> is the overall volume fraction of brushes. Q a = 
J drq a (r, s)qj,(r, s) represents the single chain partition 
function of grafted chains subject to the field w a , and 
Q/3 = J drqp(r, s)qp(r, 1 — s) is the partition function for 
solvent under the field wp. The end-segment distribu- 
tion functions qi(r, s) and q\{r, s) (or <ft(r, 1 — s)) stand 
for the probability of finding the s th segment at posi- 
tion r respectively from two ends of grafted (or free) 
chains. The qi satisfies a modified diffusion equation 
-7^ = ^f-V 2 qi — Wi(y)qi. q\ meets the same diffusion 
equation but with the right-hand side multiplied by -1. 
The interaction v is given by 



Nv 



p k B TV 



1 

V 



drlXapNtpatpp+XapNipatpp+XppNtppipp 

(4) 

In SCF theory, the fields and densities are determined 
by locating saddle points in the free energy expression 
Eq.(3). The resulting equations can be solved self- 
consistently, and here we implement the combinatorial 
screening algorithm of Drolct and Fredrickson.[57l[5^] A 
periodical boundary condition for a;-direction is applied, 
while for the z-direction, the region of z < is forbidden 
and the region of z > Z max is treated as a bath of sol- 
vent. Z max is fixed to be 80a, and the lateral simulation 
box size L x = 100a is large enough to avoid any influence 
due to the presence of colloidal particles. 

We first examine the interaction energy of two col- 
loidal particles (of radius R) deeply immersed in polymer 
brushes with particle's radii smaller than or comparable 
to polymer coils. Figure 2(a) gives the free energy F 
as a function of lateral separation d between two col- 
loidal particles, where the radii R of particles are cho- 
sen as O.LRq, 0.2.Ro, 0.25i? , 0.3i? , 0.35i? , OARq, and 
0.5i?o (from bottom to up), and the depth h is fixed at 
12.0a which is smaller than the pure brush height(around 
25.0a). i?o = aN? characterizes the natural size of poly- 
mers, and we take N — 100 here. We find that, with in- 
creasing colloidal particle's radius, the free energy curves 
change from shapes of short-range attraction and long- 
range repulsion into purely repulsive one. The energy 
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■ FIG. 2: Comparison of free energy in the present brush sys- 
tem(a) and in pure polymer melt(b), interfacial energy in the 
present brush system(c) and in pure polymer melt(d), and 
brush entropy in the present system(e) and free chain entropy 
in the pure polymer melt(f). 



barrier height is compared to thermal energy k B T, and 
the position of the peak shifts from d = 2. OA to d = 0. 
In contrast, the position of global minimum of free en- 
ergy shifts to larger values of d, indicating a transition of 
favorite colloidal aggregation to colloidal dispersion. We 
also find that there may exist a weak double-peak be- 
havior at middle particle sizes (around 0.3i?o ~ 0.35i?o), 
where a local minimum is located at d = 1.0 A and two 
peaks at d = and 2. OA. For larger particles, the free 
energy decreases with d, but at d = 2. OA, an inflexion 
point still exists. For comparison, we calculate the free 
energy curves of corresponding particles in pure poly- 
mer melts instead of polymer brush (see Fig. 2(b)). On 
the contrary, we can see that, in polymer melt, the pair 
potential between colloidal particles is monotonically at- 
tractive in spite of different sizes and the interaction is 
short-range (about 2. OA) with distance. 

In order to highlight the free energy feature shown 
in Fig. 2(a), we examine the free energy F of the col- 
loidal particles in polymer brushes which is composed 
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FIG. 3: Concentration distributions of brush chains due to the presence of two colloidal particles with size R = 0.5i?o and 
depth h = 12.0a. (a) d = 0, (b) d = 2. OA, and (c) d = 13. OA. 



of three main contributions: interfacial energy between 
colloidal particles and polymer chains, entropy contri- 
bution due to space overlap of depletion zones (deple- 
tion potential), and cntropically elastic stretching en- 
ergy of brush chains due to the compression of particles. 
The interfacial energy is the same either for the parti- 
cles immersed in polymer brushes (Fig. 2(c)) or in free 
polymer melt (Fig. 2(d)) due to the same chemical na- 
ture of grafted and free chains. The interfacial energy 
curves exhibit a turning point at d = 2. OA because the 
depletion zones surrounding the two particles will over- 
lap when the distance d is smaller than 2. OA, which will 
decrease the total interfacial area between colloidal par- 
ticles and polymer chains. Figure 2(f) shows the entropy 
curves of free polymer chains when the colloidal parti- 
cles are immersed in polymer melt. The free chain en- 
tropy change with the two-particle separation d mainly 
arises from space overlap effect of depletion zones. When 
d < 2. OA, two colloidal particles approach each other, 
and their depletion zones begin to overlap. This will in- 
crease the total volume accessible to the free polymer 
chains, leading to the increase of their entropy. Notice 
that in Fig. 2(f), slight decreases of the entropy appear 
for d > 2. OA, and the entropy curves reach their mini- 
mum at d = 3. OA. Such a behavior is due to the decrease 
of conformational entropy of polymer chains confined be- 
tween the two-particle surfaces, because when d arrives 
at 3. OA the free polymer chain begins to be permitted 
through the region between two particle' surfaces. When 
the colloidal particles are immersed in polymer brushes, 
the entropy of grafted chains will include both the contri- 
butions of space depletion effect and entropically elastic 
effect of brushes (see Fig. 2(e)). On the one hand, the 
entropy due to chain compression reaches the minimum 
at zero distance (d = 0), and increases with the distance 
d. The reason is that when the particles are far apart, 
the grafted chains will more easily fill the upper space 
through the interspace between the two-particle surfaces 
to release the chain compression. Figure 3 shows the 
brush concentration distributions due to the presence of 
two colloidal particles for three different distances d = 0, 
d = 2. OA, and d — 13. OA. The deformation of brushes 
is obviously weakened with increasing the separation of 
colloidal particles. On the other hand, the space deple- 
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FIG. 4: Free energy(a) and brush entropy(b) curves of the 
system when two particles of size R — 0.5Ro are located in 
different penetration depths. 



tion effect leads to the decrease of chain entropy with 
d for d < 2. OA where depletion zones surrounding the 
two particles overlap (see Fig. 2(f)). It is the competi- 
tion between spatial depletion effect around hard parti- 
cles and entropically elastic energy of brushes, leading 
to a possible appearance of the complex interaction for 
< d < 2. OA and the shift of minimum values of entropy 
curves from 2. OA to zero (Fig. 2(e)). For small particles, 
the space depletion effect will dominate over compression 
energy of brushes and the curve has a minimum at 2. OA. 
On the contrary, for large particles the compression effect 
of brushes will become more important and determines 
the minimum point at d = Q. It is actually interesting 
that the brushes have large entropic restoring force which 
easily controls the dispersions of colloids, in contrast to 
non-adsorbing polymer solvent. Compared Fig. 2(a) with 
Fig. 2(e), we find that the change of free energy is al- 
most attributed to the entropic effect of brushes. On the 
other hand, we see from Fig. 2(c) that the interfacial en- 
ergy mainly enhance the formation of inflexion point at 
d = 2. OA, and may lead to the onset of double peaks of 
the free energy curves for middle particle sizes. 

To understand the brush entropy effects on effective 
interactions between two colloidal particles, we now con- 
sider the influence of colloidal penetrating depth h on 
the free energy of the system for fixed radius of colloidal 
particles. Figure 4(a) and 4(b) give the free energy and 
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FIG. 5: Exponential-decay function fitting errors for ampli- 
tude A and decay length £. (a) errors for two colloidal parti- 
cles of different sizes immersed in polymer brush; (b) errors 
for two colloidal particles of size R — 0.5i?o under different 
penetrating depths h. 



brush entropy curves, respectively, for the particle size 
R = 0.5R with depths h = 20.0a, 17.0a, 15.0a, 12.0a. 
We can find that similar changes of curve shapes for 
different h values occur with the increase of d. At 
high h, entropically compression energy of brush chains 
is small compared with space depletion contributions. 
The free energy has a minimum value at d = where 
brush entropy reaches its maximum. As h is small, the 
compression contribution of brushes becomes dominant, 
and the minimum value of free energy tends to shift to 
large separation of the particles. Correspondingly, the 
minimum point of brush entropy moves from 2. OA to 
zero(Fig.4(b)). In particular, for middle depths where 
the spatial depletion effect of hard particles in polymer 
solution may compete with entropically elastic energy of 
brushes, the richer interaction potentials such as repul- 
sive barrier or double peaks may appear for free energy 
curves. 

To further examine the long-range force features due 
to entropically elastic effect of brushes, finally we de- 
fine an exponential-decay function of free energy as 
Co + A exp(-j) with increasing the two-particle separa- 
tion d. Figure 5 shows the fitting errors in percentage for 
amplitude A and decay length £ as functions of colloidal 
particle size and colloidal penetrating depth h, respec- 
tively. We find from Fig. 5 (a) that when colloidal parti- 
cles are immersed in brushes, the errors are very small 
in spite of different particle sizes. This means that the 
long-range particle-particle interaction mediated by poly- 
mer brushes has an excellent exponential-decay repulsive 
tail. Figure 5(b) provides the corresponding exponential- 
decay fitting errors in percentage for the long-range be- 
havior of free energy curves as a function of penetration 
depth h. We see that the errors, which reflect the degree 
of agreement with the exponential decay behavior, have 
a sharp transition between h — 17.0a and 20. 0a where 
the particles start to reach brush-solvent surface (see the 
inset of Fig. 5(b)). The result further addresses the fact 
that the deformation of brushes drives the occurrence of 



the exponential decay of long-range repulsions. 

In summary, we have revealed a complex interaction 
between colloidal particles mediated by polymer brushes. 
Under different colloidal sizes and penetrating depths, 
we obtain complex attractive and repulsive interaction 
behaviors as a function of colloidal particles' separa- 
tion, due to the complex interplay among entropically 
stretching energy of grafted polymers, entropic effects 
due to "depletion zone" overlap of hard particles, and 
the particles-polymer interfacial energy. The short-range 
attraction is mainly attributed to depiction effects of par- 
ticles in polymer solutions, whereas the long-range in- 
teraction between penetrating particles is found to be a 
typical result of brush entropy effect, which has a well- 
fitted exponential-decay repulsive tail. By simply using 
the grafting substrate, the present result indicates a pos- 
sibility for governing the tendency of colloidal particles 
to flocculate into multi-particle aggregates or to remain 
dispersed, and may have important implications for in- 
dustrial processes, especially for biological systems. 

Finally, we should stress that in the present case, di- 
mensionality is also an important factor. Here, the two- 
dimensional results can be effectively used to describe 
the interaction between two parallel long colloidal rods 
in three-dimensional experiment systems. However, for 
the case of colloidal spheres in three-dimensional systems, 
the deformation of grafted chains would be weakened (or 
larger colloidal spheres are needed to yield the equivalent 
deformation of grafted chains as in the two-dimensional 
case). Then, the peaks of the free energy curves(for ex- 
ample, in Fig. 2(a)) are expected to decrease under the 
same parameters, and the transition between attractive 
and repulsive interactions will finally happen at larger 
colloidal radius. Additionally, we also point out that the 
free energy difference caused by liquid and solid brushes 
will greatly depend on the grafting density of brushes. 
For low grafting densities, the chains of liquid brushes 
can easily relieve the deformation imposed by the col- 
loidal particles by moving out of the region between col- 
loids, which favors the effective attraction between col- 
loidal particles at short distances. On the contrary, for 
solid brush case where the end of chains cannot move 
away and the grafting density keeps uniform, the total 
free energy and brush entropy will depend on the relative 
lateral positions between the colloidal particles and the 
grafting end-points of chains on the substrate. For ex- 
ample, the existence of immobile grafting chain between 
the colloidal particles will enhance the entropy repulsive 
barrier to particle-particle interaction at short distances. 
This means that the resulting curves of free energy and 
brush entropy obtained for mobile and immobile grafted 
chains will be obviously different. However, at high graft- 
ing density, even the mobile grafted chains cannot move 
easily due to the incompressibility constraint, i.e., the 
grafting density almost keeps uniform on the substrate 
for mobile grafted chains. On the other hand, the brush 
height increases with increasing the grafting density, and 
thus the grafted end effect on particles becomes less im- 



6 



portant. In this case, there is no distinct difference be- This work was supported by the National Natural Sci- 
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